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(Dated: February 5, 2008) 

We present upper limits on the gravitational wave emission from 78 radio pulsars based on data 
from the third and fourth science runs of the LIGO and GEO 600 gravitational wave detectors. The 
data from both runs have been combined coherently to maximise sensitivity. For the first time 
pulsars within binary (or multiple) systems have been included in the search by taking into account 
the signal modulation due to their orbits. Our upper limits are therefore the first measured for 
56 of these pulsars. For the remaining 22, our results improve on previous upper limits by up to 
a factor of 10. For example, our tightest upper limit on the gravitational strain is 2.6 x 10^^^ for 
PSR J1603-7202, and the equatorial ellipticity of PSR J2124-3358 is less than 10"*^. Furthermore, 
our strain upper limit for the Crab pulsar is only 2.2 times greater than the fiducial spin-down limit. 

PACS numbers: 04.80.Nn, 95.55.Ym, 97.60.Gb, 07.05. Kf 



I. INTRODUCTION 

This paper details the results of a search for gravita- 
tional wave signals from known radio pulsars in data from 
the third and fourth LIGO and GEO 600 science runs (de- 
noted S3 and S4). These runs were carried out from 31*^' 
October 2003 to 9*'^ January 2004 and 22""^ February to 
23'''^ March, 2005 respectively. We have applied, and ex- 
tended, the search technique of Dupuis and Woan [l[ to 
generate upper limits on the gravitational wave ampli- 
tude from a selection of known radio pulsars, and infer 
upper limits on their equatorial ellipticities. The work is 
a natural extension our previous work given in Refs. 



A. Motivation 

To emit gravitational waves a pulsar must have some 
mass (or mass-current) asymmetry around its rotation 
axis. This can be achieved through several mechanisms 
such as elastic deformations of the solid crust or core 
or distortion of the entire star by an extremelystrong 
misaligned magnetic field (sec Sec. Ill of Ref. [j] for a 
recent review). Such mechanisms generally result in a 
triaxial neutron star which, in the quadrupole approxi- 
mation and with rotation and angular momentum axes 
aligned, would produce gravitational waves at twice the 
rotation frequency. These waves would have a charac- 
teristic strain amplitude at the Earth (assuming optimal 
orientation of the rotation axis) of 



ho 



(1.1) 



where v is the neutron star's spin frequency, Izz its prin- 
cipal moment of inertia, e = {Ixx~Iyy)/Izz^is equatorial 
ellipticity, and r its distance from Earth [5|. 

A rotating neutron star may also emit gravitational 
waves at frequencies other than 2i>. For instance, if the 
star is undergoing free precession there will be gravita- 
tional wave emission at (or close to) both i/ and 2^ Q. 



In general such a precession would modulate the time- 
of- arrival of the radio pulses. No strong evidence of such 
a modulation is seen in any of the pulsars within our 
search band, although it might go unnoticed by radio 
astronomers, cither because the modulation is small (as 
would be the case if the precession is occurring about 
an axis close to the pulsar beam axis) or because the 
period of the modulation is very long. However, this 
misalignment and precession will be quickly darnped un- 
less sustained by some mechanism (e.g. Ref. [3]), and 
even with such a mechanism calculations give strain am- 
plitudes which would probably be too low compared to 
LIGO sensitivities 0,131 • For these reasons, and for the 
reason discussed in i^IIIl we restrict our search to twice 
the rotation frequency. Of course, it cannot be ruled out 
that there are in fact other gravitational wave compo- 
nents, perhaps caused either by a stronger than expected 
precession excitation mechanism or by an event in the 
pulsar's recent past that has set it into a precessional 
motion which has not yet decayed away. A search for 
gravitational waves from the Crab pulsar at frequencies 
other than twice the rotation frequency is currently under 
way and will be presented elsewhere. 

Known pulsars provide an enticing target for gravi- 
tational wave searches as their positions and frequencies 
are generally well-known through radio or X-ray observa- 
tions. As a result the signal search covers a much smaller 
parameter space than is necessary when searching for sig- 
nals from unknown sources, giving a lower significance 
threshold. In addition, the deterministic nature of the 
waves allows a building up of the signal-to-noise ratio by 
observing coherently for a considerable time. The main 
drawback in a search for gravitational waves from the 
majority of known pulsars is that the level of emission 
is likely to be lower than can be detected with current 
detector sensitivities. 

Using existing radio measurements, and some reason- 
able assumptions, it is possible to set an upper limit on 
the gravitational wave amplitude from a pulsar based 
purely on energy conservation arguments. If one assumes 
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that the pulsar is an isolated rigid body and that the 
observed spin-down of the pulsar is due to the loss of 
rotational kinetic energy as gravitational radiation (i.e., 
dEiot/dt = A'k'^IzzVv) then the gravitational wave ampli- 
tude at the Earth (assuming optimal orientation of the 
rotation axis) would be 

Of course these assumptions may not hold, but it would 
be surprising if neutron stars radiated significantly more 
gravitational energy than this. With these uncertainties 
in mind, searches such as the one described in this paper 
place direct upper limits on gravitational wave emission 
from rotating neutron stars, and these limits are already 
approaching the regime of astrophysical interest. 



B. Previous results 

Before the advent of large-scale interfcromctric detec- 
tors, there was only a limited ability to search for gravi- 
tational waves from known pulsars. Resonant mass grav- 
itational wave detectors are only sensitive in a relatively 
narrow band around their resonant frequency and so can- 
not be used to target objects radiating outside that band. 
A specific attempt to search for gravitational waves from 
the Crab pulsar at a frequency of ^ 60 Hz was, however, 
made with a specially designed aluminium quadrupole 
antenna [1, [13] giving a la upper limit of /iq < 2x10^'^^. 
A search for gravitational waves from what was then 
the fastest millisecond pulsar, PSR J1939-I-2134, was con- 
ducted by Hough et al using a split bar detector, 
producing an upper limit of /iq < lO^^''. 

The first pulsar search using interferometer data was 
carried out with the prototype 40 m interferometer at 
Caltcch by Hercld • The search was again for gravita- 
tional waves from PSR J1939+2134, and produced upper 
limits of /lo < 3.1x10^^^ and ho < 1.5xl0~^^ for the first 
and second harmonics of the pulsar's rotation frequency. 

A much larger sample of pulsars is accessible to broad- 
band interferometers. As of the beginning of 2005 the 
Australia Telescope National Facility (ATNF) online pul- 
sar catalogue [11] listed^ 154 millisecond and young pul- 
sars, all with rotation frequencies > 25 Hz (gravitational 
wave frequency > 50 Hz) that fall within the design 
band of the LIGO and GEO 600 interferometers, and 
the search for their gravitational waves has developed 
rapidly since the start of data-taking runs in 2002. Data 
from the first science run (SI) were used to perform 
a search for gravitational waves at twice the rotation 
frequency from PSR J1939+2134 [1]. Two techniques 



The catalogue is continually updated and as such now contains 
more objects. 



were used in this search: one a frequency domain, frc- 
quentist search, and the other a time domain, Bayesian 
search which gave a 95% credible amplitude upper limit 
of 1.4x10"^^, and an eUipticity upper hmit of 2.9 xlO"'* 
assuming I^z = lO'^^kgm^. 

Analysis of data from the LIGO S2 science run set up- 
per limits on the gravitational wave amplitude from 28 
radio pulsars To do this, new radio timing data were 
obtained to ensure the pulsars' rotational phases could be 
predicted with the necessary accuracy and to check that 
none of the pulsars had glitched. These data gave strain 
upper limits as low as a few times 10~^^, and several el- 
lipticity upper limits less than 10~^. The Crab pulsar 
was also studied in this run, giving an upper limit a fac- 
tor of ~ 30 greater than the spin-down limit considered 
above. Prior to this article these were the most sensitive 
studies made. Preliminary results for the same 28 pul- 
sars using S3 data were given in Dupuis (2004) [l3| , and 
these are expanded below. 

In addition to the above, data from the LIGO S2 
run have been used to perform an all-sky (i.e., non- 
targeted) search for continuous wave signals from isolated 
sources, and a search for a signal from the neutron star 
within the binary system Sco-Xl Q . An all-sky continu- 
ous wave search using the distributed computing project 
lEinstein@homef has also been performed on S3 data [Ij] . 
These searches use the same search algorithms, are fully 
coherent and are ongoing using data from more recent 
(and therefore more sensitive) runs. Additional contin- 
uous wave searches using incoherent techniques are also 
being performed on LIGO data [H, [l^ . 

Unfortunately the pulsar population is such that most 
have spin frequencies that fall below the sensitivity band 
of current detectors. In the future, the low- frequency sen- 
sitivity of VIRGO [3 and Advanced LIGO }^ should 
allow studies of a significantly larger sample of pulsars. 



C. The signal 

Following convention, we model the observed phase 
evolution of a pulsar using a Taylor expansion about a 
fixed epoch time to- 

(t){T) = (l)o + 2^[vo{T -h) + ]^MT -tof + 

ii>o(T-to)=^ + ...}, (1.3) 

where (/jq is the initial (epoch) spin phase, i^o and its time 
derivatives are the pulsar spin frequency and spin-down 
coefficients at and T is the pulsar proper time. 

The expected signal in an interferometer from a triaxial 



^ |http : //elnsteln . phys ■ uwm . edu] 
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pulsar is 

h{t) = ^F+{t;i')ho{l + cos'^ L)cos2(j){t) + 

Fx(<; V')^ocostsin2(/)(t), (1.4) 

where (/)(i) is the phase evolution in the detector time 
t, and Fx are the detector antenna patterns for the 
plus and cross polarisations of gravitational waves, ip is 
the wave polarisation angle, and l is the angle between 
the rotation axis of the pulsar and the line-of-sight. A 
gravitational wave impinging on the interferometer will 
be modulated by Doppler, time delay and relativistic ef- 
fects caused by the motions of the Earth and other bodies 
in the solar system. Therefore we need to transform the 
'arrival time' of a wave-crest at the detector, t, to its 
arrival time at the solar system barycentre (SSB) th via 

tb = t + 6t = t+^ + AEQ+AsQ, (1.5) 

where r is the position of the detector with respect to the 
SSB, n is the unit vector pointing to the pulsar, Aeq is 
the special relativistic Einstein delay, and Asq is the 
general relativistic Shapiro delay [lOl- Although pulsars 
can be assumed to have a large velocity with respect to 
the SSB, it is conventional to ignore this Doppler term 
and set = T, as its proper motion is generally negli- 
gible (sec WI Al for cases where this assumption is not 
the case). For pulsars in binary systems, there will be 
additional time delays due to the binary orbit, discussed 

in fjnui 



H2 interferometers also both improved their duty factors 
to 80.5% and 81.4%, with longest lock stretches of almost 
a day. The typical strain sensitivities of all the interfer- 
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FIG. 1: Median strain amplitude spectral density curves for 
the LIGO and GEO 600 interferometers during the S4 run. 

omcters during S4 can be seen in Fig. [T] This shows the 
LIGO detectors reach their best sensitivities at about 
150 Hz, whilst GEO 600 achieves its best sensitivity at 
its tuned frequency of 1 kHz. 



II. INSTRUMENTAL PERFORMANCE IN S3/S4 

The S3 and S4 runs used all three LIGO interferome- 
ters (HI and H2 at the Hanford Observatory in Washing- 
ton, and LI at the Livingston Observatory in Louisiana) 
in the USA and the GEO 600 interferometer in Hannover, 
Germany. GEO 600 did not run for all of S3, but had two 
main data taking periods between which improvements 
were made to its sensitivity. All these detectors had dif- 
ferent duty factors and sensitivities. 

A. LIGO 

For S3 the HI and H2 interferometers maintained rel- 
atively high duty factors of 69.3% and 63.4% respec- 
tively. The LI interferometer was badly affected by an- 
thropogenic seismic noise sources during the day and thus 
had a duty factor of only 21.8%. 

Between S3 and S4 the LI interferometer was upgraded 
with better seismic isolation. This greatly reduced the 
amount of time the interferometer was thrown out of its 
operational state by anthropogenic noise, and allowed it 
to operate successfully during the day, with a duty factor 
of 74.5% and a longest lock stretch of 18.7 h. The HI and 



B. GEO 600 

During S3 GEO 600 was operated as a dual-recycled 
Michelson interferometer tuned to have greater sensi- 
tivity to signals around IkHz. The first period of 
GEO 600 participation in S3 was between 5*'' to U*'' 
November 2003, called S3 1, during which the detec- 
tor operated with a 95.1% duty factor. Afterwards, 
GEO 600 was taken offline to allow further commission- 
ing work aimed at improving sensitivity and stability. 
Then from 30*'^ December 2003 to 13**^ January 2004 
GEO 600 rejoined S3, called S3 II, with an improved duty 
factor of 98.7% and with more than one order of magni- 
tude improvement in peak sensitivity. During S3 there 
were five locks of longer than 24 hours and one lock longer 
than 95 hours. For more information about the perfor- 
mance of GEO 600 during S3 see Ref. [2l|. 

GEO 600 participated in S4 from 22"'^ February to 24**^ 
March 2005, with a duty factor of 96.6%. It was oper- 
ated in essentially the same optical configuration as in S3. 
With respect to S3, the sensitivity was improved more 
than an order of magnitude over a wide frequency range, 
and close to two orders or magnitude around 100 Hz. For 
more information about GEO 600 during S4 sec Ref. [l^l ■ 
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C. Data quality 

When a detector is locked on resonance and all con- 
trol loops are in their nominal running states and there 
are no on-site work activities that are known to compro- 
mise the data, then the data are said to be science mode. 
All science mode data arc not of sufficient quality to be 
analysed however, and may be flagged for exclusion. Ex- 
amples of such data quality flags are ones produced for 
epochs of excess seismic noise, and the flagging of data 
corrupted by overflows of photodiode ADCs. For this 
analysis we use all science mode data for which there is 
no corresponding data quality flag. For S3 this gives ob- 
servation times of 45.5 days for HI; 42.1 days for H2; and 
13.4 days for LI. For S4 this gives observation times of 
19.4days for HI; 22.5days for H2; and IT.ldays for LI. 

III. THE SEARCH METHOD 

Our search method involves heterodyning the data us- 
ing the phase model 4){t) to precisely unwind the phase 
evolution of the expected signal, and has been discussed 
in detail in Ref. After heterodyning, the data are 
low-pass filtered, using a ninth order Butterworth filter 
with a knee frequency of 0.5 Hz, and re-binned from the 
raw data sample rate of 16 384Hz to 1/60 Hz i.e., one 



sample per minute. The motion of the detector within 
the solar system modulates the signal and this is taken 
into account within the heterodyne by using a time de- 
lay given in Eq. (|1.5p . which transforms the signal to the 
SSB. Signals from binary pulsar systems have an extra 
modulation term in the signal, as discussed briefly below, 
and these we targeted for the first time in S3/S4. 

The search technique used here is currently only able to 
target emission at twice the pulsar's rotation frequency. 
Emission near the rotation frequency for a precessing 
star is likely to be offset from the observed pulsation fre- 
quency by some small factor dependent on unknown de- 
tails of the stellar structure . As our search technique 
requires precise knowledge of the phase evolution of the 
pulsar such an additional parameter cannot currently be 
taken into account. For the emission at twice the rotation 
frequency there is no extra parameter dependence on the 
frequency and this is what our search was designed for. 

We infer the pulsar signal parameters, denoted a = 
(/iQ, 00, cos t, V'), from their (Baycsian) posterior proba- 
bility distribution (pdf) over this parameter space, as- 
suming Gaussian noise. The data are broken up into time 
segments over which the noise can be assumed stationary 
and we analytically marginalise over the unknown noise 
fioor, giving a Student-t likelihood for the parameters for 
each segment (see Ref. [l[ for the method). Combining 
the segments gives an overall likelihood of 



M 



(3.1) 



where each is a heterodyned sample with a sample 
rate of one per minute, M is the number of segments 
into which the whole data set has been cut, rrij is the 
number of data points in the jth segment, and y/j, given 

by 

Vk = ii^+(ifc;V)/io(l+cos2 t)e''2'^"-^i^x(ifc;V')/ioCOSte'20«, 

(3.2) 

is the gravitational wave signal model evaluated at i/t, 
the time corresponding to the fc'th heterodyned sample. 
In Ref. the value of nij was fixed at 30 to give 30 min 
data segments, and data that was contiguous only on 
shorter timcscales, and which could not be fitted into 
one of these segments, was thrown out. In the analy- 
sis presented here we have allowed segment lengths to 
vary from 5 to 30 min, so we maximise the number of 30- 
minute segments whilst also allowing shorter segments at 
the end of locked stretches to contribute. The likelihood 
in Eq. p.ip assumes that the data is stationary over each 
of these 30 minute (or smaller) segments. This assump- 



tion holds well for our data. Large outliers can also be 
identified and vetoed from the data, for example those 
at the beginning of a data segment caused by the impul- 
sive ringing of the low-pass filter applied after the data 
is heterodyned. 

The prior probabilities for each of the parameters are 
taken as uniform over their respective ranges. Upper lim- 
its on /iQ a-re set by marginalising the posterior over the 
nuisance parameters and then calculating the h^'^" value 
that bounds the cumulative probability for the desired 
credible hmit of 95%: 

0.95= / ° p{K\{Bk})Aho. (3.3) 
Jo 

A. Combining data 

In the search of Ref. [H the combined data from the 
three LIGO interferometers were used to improve the sen- 
sitivity of the search. This was done by forming the joint 
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likelihood from the three independent data sets: 
p(Bfe|a)joint = p{Bk\sL)m .p{Bk\a)H2 .p{Bk\a.)Li- (3.4) 

This is valid provided the data acquisition is coherent 
between detectors and supporting evidence for this is 
presented in jVl It is of course a simple matter to ex- 
tend Eq. (j3.4p to include additional likelihood terms from 
other detectors, such as GEO 600. 

In this analysis we also combine data sets from two 
different science runs. This is appropriate because S3 
and S4 had comparable sensitivities over a large portion 
of the spectrum. Provided the data sets maintain phase 
coherence between runs, this combination can simply be 
achieved by concatenating the data sets from the two 
runs together for each detector. 

An example of the posterior pdfs for the four un- 
known pulsar parameters of PSR J0024— 7204C (each 
marginalised over the three other parameters) is shown 
in Fig. [21 The pdfs in Fig. [2] are from the joint analysis 




cosl \\f 



FIG. 2: The marginalised posterior pdfs for the four unknown 
pulsar parameters ho, c^o, cost and ^p, for PSR J0024— 7204C 
using the joint data from the three LIGO detectors over S3 
and S4. 

of the three LIGO detectors using the S3 and S4 data, 
all combined coherently. The shaded area in the ho pos- 
terior shows the area containing 95% of the probability 
as given by Eq. p.3p . In this example the posterior on 
ho is peaked at Hq = 0, though any distribution that is 
credibly close to zero is consistent with Hq = 0. Indeed 
an upper limit can formally be set even when the bulk of 
the probability is well away from zero (see the discussion 
of hardware injections in SjV|. 



B. Binary models 

Our previous known pulsar searches 0, have ex- 
cluded pulsars within binary systems, despite the ma- 
jority of pulsars within our detector band being in such 



systems. To address this, we have included an additional 
time delay to transform from the binary system barycen- 
tre (BSB) to pulsar proper time, which is a stationary 
reference frame with respect to the pulsar. The code for 
this is based on the widely used radio pulsar timing soft- 
ware TEMPO [H]. The algorithm and its testing are 
discussed more thoroughly in Ref. [2^ . 

There are five principal parameters describing a Kep- 
lerian orbit: the time of periastron, Tq, the longitude of 
periastron, ojq, the eccentricity, e, the period, Pb, and the 
projected semi- major axis, x = a sin z. These describe the 
majority of orbits very well, although to fully describe the 
orbit of some pulsars requires additional relativistic pa- 
rameters. The basic transformation and binary models 
below are summarised by Taylor and Weisberg [2Q{ and 
Lange et al. [H], and are those used in TEMPO. The 
transformation from SSB time tb to pulsar proper time 
T follows the form of Eq. (|1.5p and is 

ffc = r + Ar + Ae + As, (3.5) 

where Ar is the Roemer time delay giving the propaga- 
tion time across the binary orbit, Ae is the Einstein delay 
and gives gravitational redshift and time dilation correc- 
tions, and Ag is the Shapiro d elay and gives the general 
relativistic correction (see Ref. [20| for definitions of these 
delays). 

The majority of binary pulsars can be described by 
three orbital models: the Blandford-Teukolsky (BT) 

model, the low eccentricity (ELLl) model^ and the 

Damour-Deruelle (DD) model (see Refs. [Hl^l HI] for 
further details of these models). These different models 
make different assumptions about the system and/or are 
specialised to account for certain system features. For ex- 
ample, the ELLl model is used in cases where the eccen- 
tricity is very small, and therefore periastron is very hard 
to define, in which case the time and longitude of perias- 
tron will be highly correlated and have to be reparame- 
teriscd to the Laplace-Lagrange parameters [1^ . When a 
binary pulsar's parameters are estimated from radio ob- 
servations using TEMPO the different models are used 
accordingly. These models can be used within our search 
to calculate all the associated time delays and therefore 
correct the signal to the pulsar proper time, provided we 
have accurate model parameters for the pulsar. 

IV. PULSAR SELECTION 

The noise floor of the LIGO detectors increases rapidly 
below about 50 Hz, so pulsar targets were primarily se- 
lected on their frequency. The choice of a 50 Hz grav- 
itational wave frequency cut-off (pulsar spin frequency 
of 25 Hz) is somewhat arbitrary, but it also loosely re- 
flects the split between the population of fast (millisec- 
ond/recycled and young) pulsars and slow pulsars. 

All 154 pulsars with spin frequencies > 25 Hz were 
taken from the ATNF online pulsar catalogue [l^ (de- 
scribed in Ref. [1^). The accuracy of these parameters 
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varies for each pulsar and is dependent on the time span, 
density of observations and the noise level of the timing 
observations. Clearly it is important to ensure that pa- 
rameter uncertainties do not lead to unacceptable phase 
errors in the heterodyne. Pulsars are not perfect clocks, 
so the epoch of the parameters is also important as more 
recent measurements will better reflect the current state 
of the pulsar. Importantly, there is near-continuous mon- 
itoring of the Crab pulsar at Jodrell Bank Observatory, 
and as such its parameters are continuously updated [27| . 

Precise knowledge of the phase evolution of each target 
pulsar is vital for our analysis, and possible effects that 
may lead to a departure from the simple second-order 
Taylor expansion are discussed below. 



A. Pulsar timing 

Using TEMPO, we obtained the parameters of 75 pul- 
sars from the regular observation programs carried out 
at Jodrell Bank Observatory and the Parkes Telescope 
(see Ref. [28[ for details of the techniques used for this). 
For 37 of these the timings spanned the period of S3. 
These same model parameters were used to extrapolate 
the pulsar phases to the period of S4. The effect of pa- 
rameter uncertainties on this extrapolation is discussed 
in i jlVBi but is only important in its effect on the extrap- 
olated phase. For those pulsars observed during S3 the 
interpolation is taken to be free from significant error. 

The parameters for 16 additional pulsars (for which 
new timings were not available) were taken directly from 
the ATNF catalogue, selected using criteria described 
in the following section. The parameters of the X-ray 
pulsar PSRJ0537-6910 were taken from Ref. [H] and 
those for the Crab pulsar from the Jodrell Bank monthly 
ephemeris [27| . The remaining 61 pulsars (from the orig- 
inal list of 154) were not timed with sufficient confidence 
and were excluded from the search. This included many 
of the newly discovered pulsars (for example the 21 mil- 
lisecond pulsars in the Terzan 5 globular cluster [33|) for 
which accurate timing solutions have yet to be published. 
We therefore had a catalogue of 93 timed pulsars for our 
gravitational wave search. 



We took a 'worst-case scenario' approach by adding 
and subtracting the quoted uncertainties from the best- 
fit values of all the parameters to determine the com- 
bination which gave a maximum phase deviation, when 
propagated over the period of the run (either S3 or S4), 
from the best fit phase value calculated over the same 
time period. For example if wc assume (/'(iss) given by 
Eq. (|1.3|) (ignoring for simplicity the 00 smd v terms) 
is the best fit phase over the time span of S3, iss, the 
maximum phase uncertainty is 



A0O 



max 
1 



(^(tss) ± 27r|(i/ ± cr^)(ts3 ± crts3 



} 



(4.1) 



where the crs arc the uncertainties on the individual pa- 
rameters. Correlations between the parameters mean 
that this represents an upper limit to the maximum phase 
uncertainty, sometimes greatly overestimating its true 
value. 

There are 12 pulsars with overall phase uncertainty 
> 30° in S3, which we take as the threshold of ac- 
ceptability. A 30° phase drift could possibly give a 
factor of ~ 1 — cos 30° = 0.13 in loss of sensitiv- 
ity for a signal. Nine of these are in binary sys- 
tems (PSRsJ0024-7204H, J0407-hl607, J0437-4715, 
J142G-5625, J1518-h0205B, J1709-h2313, J1732-5049, 
J1740-5340 and J1918-G642) and in five of these Tq and 
loq contribute most to the phase uncertainty. For the 
three isolated pulsars (PSRs J0030-h0451, J0537-6910, 
and J1721— 2457) the phase error is dominated by uncer- 
tainties in frequency and/or position. 

Applying the same criterion to the time-span of S4 
we find that PSR J1730-2304 rises above the limit. For 
this pulsar its parameter uncertainties do not affect it 
for the S3 analysis as it was timed over this period, how- 
ever when extrapolating over the time of the S4 run the 
uncertainties become non-negligible. 

In total there are 13 pulsars rejected over the combined 
run. This highly conservative parameter check reduces 
our 93 candidate pulsars to 80. 



C. Timing noise 



B. Error propagation in source parameters 

The impact of parameter uncertainties on the search 
was assessed for both the S3 and S4 runs. At some level 
there are positional, frequency and frequency derivative 
uncertainties for all the target pulsars, and for pulsars 
in binary system there are also uncertainties associated 
with all the binary orbital parameters. Some of these 
uncertainties are correlated, for example the error on fre- 
quency could affect the accuracy of the first frequency 
derivative, and the binary time of periastron and longi- 
tude of periastron are also highly correlated. 



Pulsars are generally very stable rotators, but there 
are phenomena which can cause deviations in this stabil- 
ity, generically known as timing noise. The existence of 
timing noise has been clear since the early days of pul- 
sar astronomy and appears as a random walk in phase, 
frequency or frequency derivative of the pulsar about the 
regular spin-down model given in Eq. (|1.3p 13 ll . The 
strength of this effect was quantified in Ref. [31[ as an 
activity parameter A, referenced to that of the Crab pul- 
sar, and in Ref. [l^ as a stability parameter Ag. ^ is 
based on the logarithm of the ratio of the rms residual 
phase of the pulsar, after removal of the timing model, to 
that of the Crab pulsar over an approximately three- year 
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period. Ag is not based on the stochastic nature of the 
Crab pulsar's timing noise and and is defined for a fixed 
time (10* s) as 

A8=log(^i-p|x(10»s)3y (4.2) 

This assumes that the measured value of i) is dominated 
by the timing noise rather than the pulsar's intrinsic 
second spin-down derivative. Although generally true, 
this assumption is not valid for the Crab pulsar and 
PSR J0537— 6910, where a non-timing noise dominated 
i) can be measured between glitches'^. This quantity re- 
lates to the pulsar clock error caused by timing noise. 
The value of v is so small as to be unmeasurable for 
most pulsars, although an upper limit can often be de- 
fined. Arzoumanian et al. |32l | deduce, by eye, a linear 
relationship between Ag and logP of 

Ag = 6.6-l-0.61ogP, (4.3) 

where P = —v/v'^ is the period derivative. 

As defined, Ag is a somewhat imprecise indicator of the 
timing noise, not least because the time span of 10* sec- 
onds chosen by Arzoumanian et al. was simply the length 
of their data set. A preferred measure may simply be the 
magnitude and sign of P, but we shall continue to use the 
Ag parameter as our timing noise magnitude estimate for 
the current analysis. A thorough study of timing noise, 
comparing and contrasting the various measures used, 
will be given in Ref. [11] (also see Refs. [H, H^). 

There is a definite correlation between the Ag param- 
eters, spin-down rate and age. Young pulsars, like the 
Crab pulsar, generally show the most timing noise. The 
categorisation of the type of timing noise (i.e., phase, fre- 
quency or frequency derivative) in Ref. [3l| allowed them 
to ascribe different processes for each. The majority of 
pulsars studied showed frequency-type noise, possibly a 
result of random fluctuations in the star's moment of in- 
ertia. The actual mechanism behind the process is still 
unknown, with Cordes and Greenstein [35| positing and 
then ruling out several mechanisms inconsistent with ob- 
servations. 

Timing noise intrinsically linked to motions of the elec- 
tromagnetic emission source or fluctuations in the mag- 
netosphere, rather than the rotation of the pulsar, is im- 
portant in the search for gravitational waves as it may 
allow the relative phase of the electromagnetic and gravi- 
tational signals to drift. The implications of timing noise 
in this context are discussed by Jones [s^ . He gives three 
categories of timing noise, not necessarily related to the 
three types of timing noise given by Cordes and Helfand 
(3H , having different effects on any search. If all parts of 



^ These two pulsars are among the most prohfic ghtchcrs, and in 
any global fit to their parameters the value of v would most likely 
be swamped by the glitch events. 



the neutron star are strongly coupled on short timescales, 
there should be no difference between the electromag- 
netic phase and gravitational wave phase. If the tim- 
ing noise were purely a magnetospheric fluctuation then 
phase wandering caused by timing noise would not be 
seen in the gravitational wave emission. The third possi- 
bility, whereby the electromagnetic emission source wan- 
ders with respect to the mass quadrupole, could result 
from a weak exchange of angular momentum between 
the parts of the star responsible for electromagnetic and 
gravitational wave emission. Jones describes the ratio of 
the electromagnetic and gravitational timing noise phase 
residuals (A<i>) by a parameter a = A$gw/A$em, with 
the three types of timing noise described above corre- 
sponding to a = 1, and —Icm/Igw respectively, where 
the Is represent the moments of inertia of the electromag- 
netic and gravitational wave producing components. In 
principle this factor could be included as another search 
parameter. However, given the cost of including an extra 
parameter in this search, and given that it is plausible 
that that all parts of a neutron star are tightly coupled 
on the timescales of interest here, we will assume rigid 
coupling between the two components, i.e set a = 1, 
corresponding to the gravitational and electromagnetic 
signals remaining perfectly in phase. 

The Crab pulsar is regularly monitored [l^l on 
timescales that are sufficiently short to allow its timing 
noise to be effectively removed using a second heterodyne 
procedure [13] • Like the Crab pulsar, PSR J0537— 6910 is 
young, has a high glitch rate and also shows high levels of 
timing noise [2^. Unfortunately, unlike the Crab pulsar, 
we have no regular ephemeris for it that covers our data 
set, and timing irregularities are likely to be too great for 
historical data to be of use. We therefore have excluded 
PSR J0537— 6910 from the analysis. For less noisy pulsars 
we still need a method of estimating the effect of timing 
noise on phase evolution that does not rely on continu- 
ous observation. One such estimate is the Ag parameter 
given by Eq. (|4.2p . which can provide a measure of the cu- 
mulative phase error. For those pulsars with a measured 
i> we use this estimate to obtain a corresponding value of 
Ag as shown in Fig. [31 This should provide a reasonable 
estimate of the timing noise over the timespan of the pul- 
sar observation. Again we apply our criterion that cumu- 
lative phase errors of > 30° are unacceptable. In Fig. [3] 
there are four pulsars (those with the four largest Ag 
values), with measured iy, for which this is the case, and 
therefore timing noise could be a problem (having already 
noted the Crab pulsar and PSR J0537— 6910 as excep- 
tions): PSRs J1748-2446A, J1823-3021A, J1913-hl011, 
and J1952-I-3252. For pulsars with no measured i/ we 
use the approximate linear relation between the period 
derivative P and Ag given in Eq. (|4.3[) . The low P val- 
ues for these pulsars implies that timing noise will be 
negligible. 

In addition to the above, there are some pulsars in 
globular clusters for which there is no i) and for which 
P is negative {v is positive), so no value of Ag can be 
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V. HARDWARE INJECTIONS 
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FIG. 3: The values of Ag for our selection of pulsars with 
measured v. 



assigned cither through Eqs. (|4.2p or (|4.3p . For these 
pulsars the value of v (and therefore i)) must be rather 
small to have been affeeted by motions vifithin the cluster 
(discussed more in i^VI|) . so timing noise should again be 
negligible. 

For pulsars which were re-timed over the period of S3 
timing noise will be negligible (for the S3 analysis at 
least) as any timing noise, which usually has variations on 
time-scales of several months to years, will have been ab- 
sorbed in the parameter estimation. PSRs J1748— 2446A 
and J1823— 3021A were re-timed over S3 meaning that 
their S3 results will stand, although the other two will 
not. However, being conservative we will remove all four 
pulsars with large values of Ag, and PSR J0537— 6910, in 
which timing noise could be problematic, from the S4 and 
joint analysis. Note that PSR J0537— 6910 is vetoed by 
both the parameter error criterion and our timing noise 
criterion. 

This reduces our final number of well parameterised 
pulsar targets to 78 for the S3 analysis and 76 for the 
S4 and joint analyses. The 76 pulsars include 21 of 
the 28 from the previous study of Abbott et al. 3, 
and so through our selection criterion we lose the fol- 
lowing 7 previously analysed pulsars: PSRs J0030-f 0451, 
J1721-2457, J1730-2304, J1823-3021A, J1910-5959B, 
J1913-hl011, and J1952-h3252. The same selection rules 
were not applied over S2, especially of note was that no 
timing noise criterion was considered, which accounts for 
three of the pulsars we lose between the two analyses. 
Also our 30° rule was strictly applied, which the other 
four pulsars just exceeded. 

The analysis was actually performed on all 93 timed 
pulsars mentioned above, however the various parameter 
uncertainties preclude us setting upper limits on a total 
of 15 of these. 



For analysis validation purposes simulated gravita- 
tional wave signals for a variety of sources (bursts, pul- 
sars, inspirals and stochastic) have been mechanically in- 
jected into the LIGO interferometers during science runs. 
During S2 two pulsar signals were injected 0]. This was 
increased to 10 injections in the LIGO instruments for 
S3 and 12 for S4 covering a wider range of signal param- 
eters. Extracting and understanding these injections has 
been invaluable in validating the analysis. 

The hardware injection signals are produced using soft- 
ware (under LALApps [11] ) , which was largely developed 
independently of the extraction code. However, the codes 
do share the same solar system baryccntring and detec- 
tor antenna response function routines, both of which 
have been extensively checked against other sources (e.g. 
checks against TEMPO in Refs. JL|,il). 

The signals were added into each of the three LIGO 
detectors via the position control signal going to the end 
test mass in one arm. Control signals in the digital servos 
that maintain optical cavities on resonance were summed 
with fake pulsar waveforms, modulating mirror positions 
to mimic the effect of a real spinning compact object 
(i.e. differential length motions with frequency and am- 
plitude modulations appropriate for a given sky position, 
frequency and spin-down). Furthermore, as the digital 
fake waveforms have to be converted to analog coil cur- 
rents of suspended optics, the injected waveforms have to 
be divided by the transfer function of the output chain 
(predominantly the pendulum), in order to produce the 
desired differential length response of the cavity. 

The extraction of these injections is described in detail 
in Appendix]^ They show the relative phase consistency 
between the detectors over the course of a run. This 
means that a joint analysis combining the data from all 
detectors is valid. The injection plots (see Figs. [7] and [S]) 
show what we would expect our posterior plots to look 
like given a detection i.e. strongly peaked pdfs with very 
small probability at /iq = 0, as compared to those in 
Fig. [2] where ft-o peaks at zero. 



VI. RESULTS 

A. Upper limits 

Here we present 95% degree-of-belief upper limits on 
the amplitude of gravitational waves [ho) from the 78 
pulsars identified above. The value of is independent 
of any assumptions about the neutron star other than 
it is emitting gravitational waves at twice its rotation 
frequency. The results will also be presented in terms 
of the pulsars' equatorial ellipticity e, which under the 



10 



assumption of triaxiality is related to ho via Eq. (jl.ip by 

(6.1) 

To obtain an upper limit on e from that for ho we assume 
a fiducial moment of inertia value of Izz = lO'^^kgm^. 
We discuss below in Sec. IVI Bl the effect of relaxing this 
assumption. Pulsar distances are taken from the ATNF 
catalogue and are generally derived from the radio 
dispersion measures, with errors estimated to be of or- 
der 20%, although in some cases even this can be an 
underestimate. A critical review of pulsar distance mea- 
surements can be found in Ref. [3Si |. 

All upper limit results from the individual S3 and 
S4 runs along with results from the combined run, 
with and without GEO 600 included, are given in Ap- 
pendix [B] in Tables [m] and |IVl The GEO 600 data 
only provides comparable sensitivities to LIGO at fre- 
quencies greater than 1000 Hz, and are therefore only 
used in the search for PSR J1939-h2134 (at the time 
the fastest known millisecond pulsar) in S3, and addi- 
tionally PSRJ1843-1113 in S4 and the combined run. 
Inclusion of GEO 600 does not significantly change 
the joint upper limits for these pulsars. For the ma- 
jority of pulsars the lowest upper limits come from 
the combined S3/S4 data set, although for 14 pul- 
sars (PSRsJ0024-7204I, J0024-7204S, J0024-7204U, 
J0621-fl002, J1045-4509, J1757-5322 J1802-2124, 
J1804-2717, J1857-h0943, J1910-5959D, J1910-5959E, 
J1911-h0101B, J2129-5721 and J2317-hl439) the S4 re- 
sults alone provide a lower limit. The combined S3 and 
S4 run results are presented in histogram form in Fig. |4l 
Fig. [5] shows the results compared to a joint LIGO S4 
upper limit estimate curve, taken as the best sensitivity 
during S4. The joint upper limit sensitivity curve for 
the three detectors can be estimated by combining the 
detector one-sided power spectral densities (PSDs) via 

Qf f\ _ ( To\ys HI Tohs H2 Tohs LI \ /„ „n 

= 10.8^^, 

where Sh{f) is the PSD and Tobs is each detector's live 
time (using the associated duty factor of each interferom- 
eter during the run). The factor of 10.8 is given in Ref. [l| 
and was calculated through simulations with Gaussian 
noise^. 

The results are also compared to the upper limit de- 
duced from the observed spin-down via Eq. (|1.2p making 



the assumption that all rotational energy is lost through 
gravitational wave emission. The spin-down limit is seen 
as a natural crossing point after which gravitational- 
wave data, including upper limits, have a likely bear- 
ing on the nature of the neutron star. The spin-down 
upper limit will obviously depend on v. This value how- 
ever, can be masked by radial and transverse motions 
of the object (see Ref. [i^l for discussion of these ef- 
fects). The Shklovskii effect [4l|, in which the pulsar 
has a large transverse velocity will cause an apparent 
rate of change in the pulsar's period of 

Ps = -P- (6.3) 

rc 

Its 1/r dependence makes this effect mo re p rominent for 
nearby pulsars. In the ATNF catalogue [IjI values of the 
intrinsic period derivative Pint = P — Ps can be obtained 
where this effect has been corrected for. This provides a 
measure of intrinsic (rather than apparent) spin-down"" 
and, when available, is used in the spin-down ratio re- 
sults. 

The observed value of Pobs will also differ from its in- 
trinsic value. Pint, if the pulsar is accelerating - a likely 
scenario in the gravitational field of a globular cluster 
[40| . Any line-of-sight component to the acceleration, 
a|[, will give an observed value of 

Pobs = Pint + ^P (6.4) 

c 

where P is the spin period [13] ■ These effects can cause 
pulsars to have apparent spin-ups (seen in quite a large 
number of globular cluster pulsars) although are only 
strong enough to greatly affect pulsars with intrinsically 
small period derivatives. There are still many globular 
clusters for which the radial accelerations have not been 
measured, therefore no firm spin-down upper limit can be 
set, making the direct gravitational wave results a unique 
limit. 

Highlights of the combined S3/S4 results include the 
tightest strain upper limit set on a pulsar of /iq^^° = 2.6x 
10-25 for PSR J1603-7207, the smallest ellipticity at e = 
7.1x10"'^ for PSRJ2124-3358, and the closest upper 
limit to the spin-down limit at a ratio of 2.2 for the Crab 
pulsar (PSRJ0534-h2200). 



B. Dependence on the moment of inertia 

The pulsar ellipticity results detailed above assume a 
moment of inertia of lO'^^kgm'^, which is the standard 



* In Ref. a similar plot to Fig. [5] is shown for the S2 data using 
a factor of 11.4 in the relation between the upper limit and PSD. 
This definition comes from using the .7^-statistic search method 
and setting a 1% false alarm rate and 10% false dismissal rate 
for signals given the underlying detector PSD Q| 



^ Note that the heterodyne procedure still needs to make use of the 
measured spin-down rather than the intrinsic spin-down as these 
Doppler eff'ects will have the same effect on the gravitational 
waves. 
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FIG. 4: Histograms of the log of amplitude, ellipticity, and ratio of spin-down to gravitational wave upper limits for the 
combined LIGO S3 and S4 run. 
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FIG. 5: The combined S3 and S4 upper limit results on the 
amplitude of gravitational waves for 76 pulsars using LIGO 
data compared to the joint sensitivity curve. 

fiducial number used in the literature. However, modern 
theoretically computed equations of state (EOS) gener- 
ally predict somewhat larger moments of inertia for stars 
more massive than 1 Mq , a group which includes all neu- 
tron stars with measured masses (see Ref. [1^). There- 
fore the dependence on moment of inertia should be con- 
sidered. 

Bejger, Bulik, and Hacnscl H^l give an overview of 
the theoretical expectations for the moment of inertia. 
Their Fig. 2 plots the moment of inertia vs. mass for 
several theoretically predicted types of EOS. The maxi- 
mum moment of inertia they find (after varying the mass 
of the star) is 2.3 times the fiducial value, with stars of 
IAMq having moments of inertia 1.2-2.0 except for one 
outlying type of EOS. Typically the maximum moment 
of inertia occurs for a neutron star mass of 1.7 A/q or 



more. Recently masses greater than 1.6-1.7Mq with 
95% confidence have been measured [13, for some 
systems, making this reasonable to consider. More re- 
cently Lackey |46| found the highest moment of inertia 
to be 3.3x10^8 kgm^ for EOS G4 of Lackey, Nayyar, and 
Owen [13]. This is a relativistic mean-field EOS simi- 
lar to the GNH family considered by Bejger, Bulik, and 
Haensel [31 but contains no exotic phases of matter such 
as hyperons or quarks. Consequently, wc consider the 
range of theoretically predicted moments of inertia to be 
approximately 1-3 x 10'^* kg m^ . 

There have been recent attempts to infer neutron 
star moments of inertia from observations. Bejger and 
Haensel [1^, HI] derived a value for the Crab pulsar's 
moment of inertia by equating the spin-down power to 
the observed electromagnetic luminosity and inferred ac- 
celeration of the nebula. However, this (extremely high) 
value is dominated by the assumptions about the highly 
uncertain mass and mass distribution of the nebula as 
well as the relativistic wind from the pulsar, and thus 
cannot yet be considered to give a reliable value. The 
double pulsar system J0737-3039 shows great promise for 
tighter measurements of the moment of inert ia ( and con- 
straints on the EOS) in the near future [4l,[53,[5l|,|53. 
However, for the moment we are left with the theoretical 
range quoted above. 

As suggested in Ref. [s^, instead of using Eq. (|6.ip 
to set a limit on s assuming a value of Izz, one can use 
it to set a limit on the neutron star quadrupole moment 
~ Izz£ without relying on any assumption about Izz ■ The 
limit on the quadrupole moment can then be used to help 
define an exclusion region in the I-e plane. This exclusion 
region allows one to read off an upper limit on e as a 
function of the EOS-dependent moment of inertia. The 
spin-down can also be used to provide exclusion regions 
via the relation 

5 I z-' [ 1 



12 



Theoretical contributions to the exclusion regions come 
from predictions of the maximum moment of inertia and 
ellipticity. In terms of the exclusion region, our observa- 
tional upper limits on ho are far from contributing except 
for the Crab pulsar, to which we now turn. 



C. The Crab pulsar - PSR J0534+2200 

Of the known radio pulsars, the Crab pulsar has of- 
ten been considered one of the most promising sources 
of gravitational waves. This is due to its youth and 
large spin-down rate, leading to a relatively large spin- 
down upper limit several orders of magnitude higher than 
for most other pulsars. The high rate of glitching in 
the pulsar also provides possible evidence of asymmetry. 
One glitch model favoured for the Crab pulsar involves 
a change in the pulsar ellipticity, and breaking of the 
crust, as the star settles to its new equilibrium state as it 
spins down [i^. In the 1970s, estimates of gravitational 
wave strains were spurred on by the experimenters pro- 
ducing novel technologies which allowed the possibility 
of probing these low strains, with Zimmermann [s^ l pro- 
ducing estimates of gravitational wave strains from the 
Crab pulsar ranging from Zip ~ 2 x 10~^^ — 10"^^. 

The first searches for gravitational waves from the 
Crab pulsar were carried out using specially designed 
resonant bar detectors, with frequencies of around 60 Hz 
0. The most recent result using such a bar was from 
1993 and gave a la upper limit of /iq < 2 x 10~^^ p^ . 
This upper limit was passed in the LIGO S2 run, which 
gave hl^°''" 4.1x10-23 [|. Using Eq. JL!]), and taking 
= lO^^kgm^ and r = 2kpc, gives a spin-down up- 
per limit for the Crab pulsar of ho < 1.4 xlO-^*, about 
a factor of 30 below the S2 observational upper limit. 
However, the S2 limit on the Crab was at the time the 
closest approach to the spin-down limit obtained for any 
pulsar. 

Our new results for the Crab pulsar (and the other 77 
targets) are shown in Table HITl The results improve by 
up to an order of magnitude over those from the S2 run, 
and the majority of this improvement was between the 
S2 and S3 runs. The results for the Crab pulsar over the 
S2, S3 and S4 runs are plotted on the I-e plane in Fig.[Sl 
The diagonal lines in Fig. [S] mark the lower boundaries of 
exclusion regions on this plane using our upper limits ob- 
tained for the different runs. The dashed black diagonal 
line marks the lower boundary of the upper limit from 
spin-down as given in Eq. (|6.5p . The dashed horizontal 
black lines give lower and upper bounds on the moment 
of inertia of l-3xl03^kgm2; as given by our arguments 
in WIBI It can be seen that our experimental results 
currently only beat the spin-down limit for moments of 
inertia at values greater than almost double the maxi- 
mum of our theoretical range. However, over this range 
the ratio of the gravitational to spin-down upper limit 
ranges from 2.2 at the lowest value to only 1.3 at the 
largest value. 















spir 


-downf 




^^SjT"^^ 
















lU 


10 ^ 

ellipticity e 


10"^ 


IC 



FIG. 6: The moment of inertia-ellipticity plane for the Crab 
pulsar over the S2, S3 and S4 runs. The areas to the right 
of the diagonal lines are the experimentally excluded regions. 
The horizontal lines represent theoretical upper and lower lim- 
its on the moment of inertia as mentioned in ijVI Bl Theoret- 
ical upper limits on the ellipticity are much more uncertain, 
the highest being a few times 10^^. 

The spin-down limit in fact over-estimates the 
strongest possible signal because we know that much of 
the spin-down energy of the Crab goes into powering the 
nebula through electromagnetic radiation and relativistic 
particle winds. Thus it is interesting to ask how far we 
would need to beat the spin-down limit to have a chance 
of detecting a signal allowing for what is known about 
the non gravitational wave spin-down. Palomba (ssj uses 
the observed braking index 2.51 of the Crab pulsar with 
a simple model of spin-down through gravitational radia- 
tion (braking index 5) combined with some other mecha- 
nism (braking index a free parameter) to place an upper 
limit of about e < 3 x lO"''. This is about 2.5 times 
lower than the spin-down limit and 5.5 times lower than 
our result (for I^^, = lO^^kgm^). 

The Crab pulsar experienced two glitches between S3 
and S4, a large glitch on 6"^ September 2004 and a 
smaller ghtch on 22"^^ November 2004 [13 . The effect 
of glitches on the relative phase between the electromag- 
netic pulse and any possible gravitational wave signal 
is unknown, so there is uncertainty whether the (phase- 
coherent) combined S3/S4 result is valid. The combined 
result stands, but the reader should be aware that it in- 
cludes the assumption of trans-glitch phase coherence. 

VII. ASTROPHYSICAL INTERPRETATION 

We have produced new, tight, upper limits on grav- 
itational wave signal strength from a large selection of 
known pulsars and for the Crab pulsar we are very near 
the fiducial limit set by spin-down arguments. 

It can be seen from Table Hill and Fig. [5] that, for the 
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majority of pulsars, the gravitational wave detector up- 
per limits are at least one hundred times above those 
from the spin-down argument, so is there anything that 
we can take from the results in terms of astrophysics? 

First we should note that spin-down limits on grav- 
itational wave luminosity are plausible, but model de- 
pendent. They assume a model for the structure of the 
neutron star (for instance, that it is not accreting and 
is rigidly rotating, in addition to assumptions about its 
equation of state) and they take dispersion measure dis- 
tance as a consistently good measure of true distance. 
There is some considerable uncertainty associated with 
all of these assumptions. In contrast our observations set 
direct limits on a source's gravitational wave strain. 

Secondly, for globular cluster pulsars the spin-down 
measured from radio timing observations is a combina- 
tion of the spin-down intrinsic to the pulsar and acceler- 
ation along the line-of-sight a|| in the cluster's gravita- 
tional potential (see Eq. (|6.4p ). In general, the magnitude 
and sign of the acceleration is unknown but the intrinsic 
Pint > of millisecond pulsars is usually small and often 
smaller than the extrinsic contribution. Only if Pobs < 
one can be sure that a\\ < 0. Therefore, the limits derived 
from our gravitational wave observations provide the only 
direct limits on Pint which are independent from biasing 
kinematic effects. These can be combined with the ob- 
served spin-down to provide a limit on the acceleration 
in the cluster, i.e. a\\ > c(Pobs - Pg™")/P. 

Finally, it is interesting to note that our cUipticity lim- 
its are well into the range permitted by some models 
of strange quark stars or hybrid stars (e ~ a few times 



10 ^ — 10 ^) and are reaching into the range permitted 
by more conventional neutron star EOSs (e ^ a few times 

lO"'^) m. 

Currently the fifth LSC science run (S5) is underway, 
and this promises to beat the Crab pulsar spin-down limit 
within a few months of its start. For many other pulsars 
we should be able to reach amplitude upper limits of 
< 1 X 10~^^ and ellipticities of ~ 1 x 10~^. 
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more recent versions of the detector calibrations, the re- 
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pulsars signals were injected, with a further one added 
at the end of the run to be in coincidence with a single 
injection into GEO 600 [s^l- The majority of injection 
parameters were decided upon randomly, although pul- 
sar frequencies were chosen to avoid major instrumental 
or calibration lines, and amplitudes were dependent on 
the frequency. The injections were split into two groups 
of five, where values of were calculated to give two each 
with signal-to-noise ratios of approximately 3, 9, 27, 81 
and 243. The parameter values are shown in TableHl The 
10 initial signals were injected into the LIGO detectors 
for approximately the first half of the run, then turned 
off for two weeks, to ensure data were present that was 
not artificially contaminated, and then turned back on 
with the two loudest signals removed. The simultaneous 
injection with GEO 600 was switched on near the end of 
the run. 

These signals were extracted from the data using the 
analysis techniques described in i ^IIII and Ref. [l| . The 
two most important parameters for checking that the 
calibration of the instruments was correct were the am- 
plitude and initial phase, so in the Bayesian parameter 
estimation procedure the t and ^ parameters were held 
fixed at their known values. This was done because the 
correlations between ho and cost and and ip, respec- 
tively could lead to the marginalised posterior probability 
distributions (pdfs) for each parameter being distorted or 
spread out (see Ref. [31 for examples of this). The ex- 
tracted pdfs of ho and 0o for each of the injections, after 
corrections described below in 3[ can be seen in Fig. [T] 

For the vast majority of signals the extracted pdfs over- 
lap with the injected value. For the strongest injections 
with the largest signal-to-noise ratios the pdfs are rather 
narrow, and any uncertainties in the calibration become 
evident, with a maximum offset of the order of 10-15%. 
The far wider pdfs associated with the LI signal injec- 
tions reflect the lower LI sensitivity and lower duty factor 
compared with the HI and H2 detectors. It can be seen 
that the injected phases for each detector agree with each 
other to within a few degrees and are within the uncer- 
tainty of the method. This provides some evidence that 
there is phase coherence between the detectors and that a 
joint analysis, combining the data from all the detectors, 
is possible. 

Two main discrepancies have been identified as oper- 
ational mistakes made during the injection procedure: 
Pulsar? was injected into 112 with a much lower ampli- 
tude than intended, and remained undetected, therefore 
no joint analysis was performed, and PulsarO was in- 
jected into HI with an amplitude 1.6 times larger than 
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TABLE I: The parameter values for the pulsar hardware injections in S3 and S4. 



Pulsar 


a (rads) 


S (rads) 


i/gw (Hz 





1.25 


-0.98 


265.5 


1 


0.65 


-0.51 


849.1 


2 


3.76 


0.06 


575.2 


3 


3.11 


-0.58 


108.9 


4 


4.89 


-0.21 


1430.2 


5 


5.28 


-1.46 


52.8 


6 


6.26 


-1.14 


148.7 


7 


3.90 


-0.36 


1221.0 


8 


6.13 


-0.58 


194.3 


9 


3.47 


1.32 


763.8 


Geo 


0.78 


-0.62 


1125.6 



i>gw (Hz/s) 
-4.15 X 10"^^ 
-3.00X 10"^° 
-1.37X 10"^^ 
-1.46X 10"^'' 
-2.54X 10"** 
-4.03X 10"^** 
-6.73x10"" 
-1.12x10"" 
-8.65x10"" 
-1.45 X 10"^'' 
-2.87x10"^^ 



ho (S3) ho (S4) (po (rads) l (rads) i/f (rads) 



4.93x10"^^ 
4.24X 10"^* 
8.04 X 10"^* 
3.26x10"^^ 
4.56x10"^^ 
9.70x10"^'' 
2.77X lO"^"* 
1.32x10"^^ 
3.18x10"^^ 
8.13x10"^'' 



9.38X lO"^'^ 
8.49 x 10"^* 
1.56 X 10"^^ 
6.16X 10~^^ 

1.01 X 10"^"^ 
1.83 X 10 
5.24 X 10 
2.81 X 10"^^ 

6.02 X 10"^^ 
1.61 X 10"^^ 
7.5x10"^^ 



-23 
-24 



2.66 
1.28 
4.03 
5.53 
4.83 
2.23 
0.97 
5.24 
5.89 
1.01 
1.99 



0.65 
1.09 
2.76 
1.65 
1.29 
1.09 
1.73 
0.71 
1.50 
2.23 
0.84 



0.77 
0.36 
-0.22 
0.44 
-0.65 
-0.36 
0.47 
0.51 
0.17 
-0.01 
0.37 



intended. 

The injection of the signal into GEO 600 is described in 
Ref. [S^l , and its analysis is described in Ref . . It was 
found that the injection performed during S3 was badly 
contaminated and could not be used. However, a sub- 
sequent injection performed shortly after S3 has verified 
that the signal parameters were correctly injected and 
extracted, validating the injection hardware and analysis 
software. 



2. S4 injections 

The 10 injections used in S3 were used again for S4 
to create artificial signals in the LIGO interferometers. 
However, their amplitudes were adjusted to give approx- 
imately the same signal-to- noise ratios as seen in S3, tak- 
ing into account the better sensitivity during the S4 run. 
For all except PulsarQ the ho values were reduced by 
a ~ half, with Pulsar9 being so strong that its ampli- 
tude was reduced by a factor of ^ 20. These signals 
were injected for the second half of the run from 8*'' 
March 2005 onwards. The updated ho values are shown 
in Table HI There were also an additional two signals 
(PulsarIO and 11), simulated to be from pulsars in bi- 
nary systems, injected for the last day of the run. The 
binary pulsar injections allowed the testing of the binary 
timing code described in ijlHBI as the injection code and 
extraction code were written independently. The binary 
injection signal parameters for PulsarIO and PulsarII 
were taken from Pulsar3 and 8 respectively, with the 
frequencies changed, and amplitudes increased to make 
sure they were visible over the short injection time. The 
frequency, amplitude and binary system parameters are 
shown in Table |TT1 The binary system parameters were 
chosen to have one in a relatively eccentric orbit and one 
in a circular orbit. We chose fairly short periods, so that 
they would have completed or nearly completed at least 
one full orbit during the injection. The To values are 
given in the pulsar rest frame. For the recovery of the bi- 
nary system injections the BT model was used, although 
as no relativistic parameters were included any of the 



models could have been used. 

The extracted amplitude and phase pdfs, after correc- 
tions described below in 3[ are shown in Fig. [5] The 
observed phase consistency between the detectors, means 
that joint hkelihoods, using all three detectors, can be 
calculated. In general the values of ho are well matched 
with the injection values. It can again be seen that for 
the strongest signals the narrow pdfs are offset from the 
injected value in ho reflecting the calibration uncertain- 
ties of 5-10%. 

The binary pulsar injections show matches to their in- 
jected values. This is a good confirmation that the binary 
timing code can track the phase well and has no signifi- 
cant errors. 



3. Calibration Issues 

A brief note should be made of the effect of calibra- 
tions on the above extracted pulsar hardware injections. 
The injections were for the most part analysed using ex- 
actly the same pipeline as applied to the general known 
pulsar analysis. However due to the nature of the hard- 
ware injections some additional post-processing of the 
results has had to be applied. To calculate the ampli- 
tude and phase of the injections, when applying forces 
to the interferometer end test masses, a reference cali- 
bration must be used. These reference calibrations are 
different for each interferometer. For both S3 and S4 
these reference calibrations differed by small, but not in- 
significant amounts, from the final calibration used when 
extracting the signals, meaning that upon extracting the 
signals the amplitude and phase appear offset from the 
injected values. As the difference between the reference 
and final calibrations are different for each interferometer 
there will also be slight offsets between the extracted pa- 
rameters between detectors. The extracted signals from 
each interferometer therefore have had to be adjusted to 
reflect these differences, determined independently of the 
hardware injections, and correct them so as to give the 
same input signal. This allows the combined joint upper 
limits to be produced. 
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TABLE II: The parameter values for the S4 binary pulsar hardware injections. 



Pulsar 


i/gw (Hz) 


ho 


To (MJD) 


Pb (days) 


e 


ujo (dcgs) 


a sin z (sees) 


10 


250.6 


1.30x10"^^ 


51749.71156482 


1.35405939 


0.0 


0.0 


1.65284 


11 


188.0 


5.21 X 10"^^ 


52812.92041176 


0.31963390 


0.180567 


322.571 


2.7564 




FIG. 7: The pdfs of ho and (j}o for 10 isolated pulsar injections into the LIGO detectors during S3. The anomaly seen in 
Pulsar7 is discussed in the text. 
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APPENDIX B: TABLES OF UPPER LIMIT RESULTS 



TABLE IIL Pulsar upper limits using LIGO data from the S3 and S4 runs. The approximate pulsar spin frequencies and spin-down rates 
globular cluster pulsars for which no spin-down upper limit could be set. The values marked with a f represent pulsars for which the spin-down 
for the Shklovskii effect. The ratio column gives the ratio of our experimental upper limits to the spin-down upper limits. 



are given. A * denotes 
limit has been corrected 
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-22 


24 


-22 


70 


-23 


24 


-2 


23 


13340 


-23 


40 


-23 


13 


-23 


36 


-23 


68 


-2 


66 


4920 


-23 


57 


-23 


13 


-23 


37 


-23 


69 


-2 


68 


4768 


J1823-3021A 


367.65 


-1 


14 X 


10" 


13 


-23 


74 


-23 


46 


-23 


36 


-23 


77 


-4 


02 


671 














































J1824-2452 


327.41 


-1 


74 X 


10" 


13 


-23 


93 


-23 


41 


-23 


05 


-23 


92 


-4 


88 


321 + 


-23 


87 


-23 


77 


-23 


68 


-24 


09 


-5 


06 


214+ 


-24 


03 


-23 


87 


-23 


70 


-24 


19 


-5 


16 


171 + 


J1843-1113 


541.81 


-2 


82 X 


10" 


15 


-23 


24 


-23 


14 


-22 


77 


-23 


30 


-5 


09 


5429 


-23 


59 


-23 


57 


-23 


43 


-23 


85 


-5 


65 


1508 


-23 


61 


-23 


61 


-23 


42 


-23 


89 


-5 


69 


1370 


J1857+0943 


186.49 


-6 


20 X 


10" 


16 


-23 


72 


-23 


50 


-23 


41 


-23 


71 


-4 


92 


1223+ 


-23 


97 


-23 


80 


-24 


04 


-24 


31 


-5 


52 


309+ 


-24 


00 


-23 


82 


-24 


07 


-24 


30 


-5 


51 


313 + 


J1905+0400 


264.24 


-3 


39 X 


10" 


16 


-23 


89 


-23 


22 


-23 


27 


-23 


83 


-5 


17 


2185 


-23 


88 


-23 


72 


-23 


81 


-23 


88 


-5 


22 


1937 


-24 


05 


-23 


84 


-23 


84 


-24 


02 


-5 


36 


1410 


J1909-3744 


339.32 


-1 


61 X 


10" 


15 


-23 


82 


-23 


41 


-22 


97 


-23 


92 


-5 


55 


1725+ 


-23 


84 


-23 


89 


-23 


67 


-24 


00 


-5 


63 


1437 + 


-23 


88 


-23 


97 


-23 


72 


-24 


13 


-5 


76 


1060 + 


J1910-5959A 


306.17 


-2 


88 X 


10" 


16 


-23 


93 


-23 


35 


-23 


10 


-24 


01 


-5 


01 


4950 


-23 


79 


-23 


59 


-23 


76 


-24 


07 


-5 


07 


4330 


-24 


01 


-23 


80 


-23 


78 


-24 


21 


-5 


20 


3178 


J1910-5959B 


119.65 


+ 1 


14 X 


10" 


14 


-24 


28 


-23 


04 


-23 


36 


-24 


29 


-4 


47 


* 


-23 


61 


-24 


02 


-24 


06 


-24 


20 


-4 


38 


* 


-24 


30 


-24 


02 


-24 


07 


-24 


39 


-4 


57 


* 



S3 and S4 



log £ I ratio 



TABLE [ml cont. 



S3 



S4 



S3 and S4 



log £ ratio 



log 



log e ratio 



Pulsar 


ly (Hz) 


l> 


(Hzs-1 


) 


HI 


H2 


LI 


Joint 


HI 


H2 


LI 


Joint 


HI 


H2 


LI 


.Joint 


J1910-5959C 


189.49 


-7 


90 X 


10" 


17 


-24 


03 


-23 


60 


-23 


37 


-24 


04 


-4 


62 


7018 


-24 


18 


-24 


07 


-23 


97 


-24 


23 


-4 


81 


4470 


-24 


32 


-24 


10 


-24 


00 


-24 


39 


-4 


97 


3107 


J1910-5959D 


110.68 


-1 


18 X 


10~ 


14 


-24 


20 


-23 


35 


-23 


26 


-24 


12 


-4 


23 


369 


-24 


40 


-24 


05 


-24 


17 


-24 


48 


-4 


60 


158 


-24 


41 


-24 


07 


-24 


12 


-24 


42 


-4 


53 


182 


J1910-5959E 


218.73 


+2 


09 X 


10" 


14 


-23 


98 


-23 


71 


-23 


30 


-24 


05 


-4 


75 


+ 


-24 


09 


-23 


96 


-23 


99 


-24 


35 


-5 


05 


+ 


-24 


15 


-24 


01 


-23 


99 


-24 


34 


-5 


04 


* 


J1911+0101A 


276.36 


+5 


03 X 


10" 


16 


-23 


56 


-23 


38 


-23 


05 


-23 


79 


-4 


43 


+ 


-23 


89 


-23 


83 


-23 


83 


-24 


12 


-4 


76 


+ 


-24 


07 


-23 


84 


-23 


81 


-24 


22 


-4 


86 


* 


J1911+0101B 


185.72 


+6 


90 X 


10" 


17 


-23 


70 


-23 


43 


-23 


30 


-23 


78 


-4 


07 


+ 


-24 


05 


-24 


01 


-24 


05 


-24 


24 


-4 


53 


+ 


-23 


96 


-24 


08 


-23 


99 


-24 


13 


-4 


43 


* 


J1911-1114 


275.81 


-1 


08 X 


10~ 


15 


-23 


85 


-23 


51 


-23 


00 


-23 


84 


-5 


14 


2188+ 


-23 


99 


-23 


80 


-23 


71 


-24 


07 


-5 


37 


1288+ 


-24 


04 


-23 


86 


-23 


70 


-24 


12 


-5 


43 


1139+ 


J1939+2134 


641.93 


-4 


33 X 


10~ 


14 


-23 


17 


-23 


22 


-22 


51 


-23 


37 


-5 


05 


2323+ 


-23 


41 


-23 


44 


-23 


33 


-23 


68 


-5 


36 


1146+ 


-23 


46 


-23 


49 


-23 


31 


-23 


78 


-5 


47 


899 + 


J1955+2908 


163.05 


-7 


91 X 


10" 


16 


-23 


82 


-23 


61 


-23 


35 


-23 


88 


-4 


20 


4088+ 


-24 


22 


-23 


94 


-24 


02 


-24 


32 


-4 


64 


1487 + 


-24 


26 


-23 


96 


-23 


99 


-24 


40 


-4 


72 


1223 + 


J1959+2048 


622.12 


-6 


52 X 


10" 


15 


-23 


30 


-23 


04 


-22 


65 


-23 


35 


-5 


38 


3220+ 


-23 


66 


-23 


37 


-23 


54 


-23 


79 


-5 


82 


1163 + 


-23 


68 


-23 


36 


-23 


54 


-23 


82 


-5 


85 


1080 + 


J2019+2425 


254.16 


-4 


54 X 


10~ 


16 


-23 


79 


-23 


49 


-23 


23 


-23 


93 


-5 


41 


1628+ 


-23 


93 


-23 


89 


-23 


79 


-24 


16 


-5 


64 


970+ 


-23 


93 


-23 


91 


-23 


75 


-24 


17 


-5 


65 


950 + 


J2051-0827 


221.80 


-6 


27 X 


10" 


16 


-23 


97 


-23 


49 


-23 


16 


-24 


05 


-5 


26 


852 + 


-23 


83 


-23 


79 


-23 


77 


-24 


02 


-5 


23 


917+ 


-24 


00 


-23 


83 


-23 


77 


-24 


11 


-5 


32 


740+ 


J2124-3358 


202.79 


-8 


45 X 


10" 


16 


-24 


08 


-23 


56 


-23 


35 


-24 


06 


-5 


90 


165 + 


-23 


91 


-23 


91 


-24 


00 


-24 


13 


-5 


98 


139+ 


-24 


20 


-24 


03 


-23 


99 


-24 


31 


-6 


15 


93 + 


J2129-5721 


268.36 


-1 


49 X 


10" 


15 


-23 


82 


-23 


45 


-23 


22 


-23 


89 


-4 


96 


1809+ 


-24 


02 


-23 


95 


-23 


90 


-24 


19 


-5 


27 


892+ 


-24 


01 


-23 


91 


-23 


89 


-24 


18 


-5 


25 


926 + 


J2140-2310A 


90.75 


+4 


27 X 


10~ 


16 


-24 


25 


-23 


48 


-23 


52 


-24 


23 


-3 


82 


* 


-24 


31 


-24 


00 


-24 


06 


-24 


45 


-4 


04 


+ 


-24 


45 


-23 


99 


-24 


10 


-24 


55 


-4 


14 


* 


J2145-0750 


62.30 


-1 


15 X 


10~ 


16 


-24 


14 


-23 


18 


-23 


53 


-24 


18 


-4 


69 


326 + 


-24 


23 


-23 


76 


-24 


24 


-24 


45 


-4 


97 


173+ 


-24 


31 


-23 


78 


-24 


26 


-24 


47 


-4 


98 


167 + 


J2229+2643 


335.82 


-1 


65 X 


10" 


16 


-23 


58 


-23 


31 


-23 


08 


-23 


64 


-5 


17 


5747 


-23 


77 


-23 


51 


-23 


72 


-23 


89 


-5 


41 


3253 


-23 


89 


-23 


56 


-23 


71 


-24 


13 


-5 


65 


1885 


J2317+1439 


290.25 


-2 


04 X 


10" 


16 


-23 


41 


-23 


40 


-23 


17 


-23 


55 


-4 


82 


9996+ 


-23 


89 


-23 


91 


-23 


79 


-24 


17 


-5 


44 


2406 + 


-23 


86 


-23 


90 


-23 


80 


-24 


15 


-5 


43 


2500 + 


J2322+2057 


207.97 


-4 


20 X 


10" 


16 


-23 


86 


-23 


48 


-23 


17 


-24 


07 


-5 


44 


900 + 


-24 


01 


-23 


91 


-23 


85 


-24 


19 


-5 


56 


673+ 


-24 


07 


-23 


96 


-23 


96 


-24 


26 


-5 


63 


578 + 



TABLE IV: Upper limits including GEO 600 for the two fastest pulsars in the analysis. 





S3 


S4 


S3 and S4 




log/iQ^'^ 1 logs 1 ratio 


log/iQ^'^ 1 logs 1 ratio 


log/iQ''^'* 1 log £ 1 ratio 


Pulsar 


f (Hz) 


l> (Hzs"l) 


HI 


H2 


LI 


GEO 


Joint 


HI 


H2 


LI 


GEO 


Joint 


HI 


H2 


LI 


GEO 


Joint 


J1843-1113 
J1939+2134 


541.81 
641.93 


-2.82X10"'^ 
-4.33X10"" 


-23.24 
-23.17 


-23.14 
-23.22 


-22.77 
-22.51 


+ 

-22.23 


-23.30 
-23.37 


-5.09 
-5.05 


5429 
2323+ 


-23.59 
-23.41 


-23.57 
-23.44 


-23.43 
-23.33 


-22.45 
-22.68 


-23.85 
-23.66 


-5.65 
-5.34 


1508 
1189 + 


-23.61 
-23.46 


-23.61 
-23.49 


-23.42 
-23.31 


-22.45 
-22.71 


-23.90 
-23.77 


-5.67 
-5.46 


1348 
914+ 



